Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 



(MN MfeK style file vl.4) 



Non-Linear Effects on the Angular Correlation Function 



Ryan Scranton 1,2 and Scott Dodelson 1 ' 2 

1 NASA/Fermilab Astrophysics Center, P.O. Box 500, Batavia, IL 60510 USA 

2 Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637 USA 



1 February 2008 



ABSTRACT 



o 
o 
o 

(N 



> 

m 
O 
m 

o 
o 
o 



Oh' 
I 1 

o 

CO 

c3 



X 



Extracting the three dimensional power spectrum from the 2D distribution of 
galaxies has become a standard tool of cosmology. This extraction requires some as- 
sumptions about the scaling of the power spectrum with redshift; all treatments to 
date assume a simple power law scaling. In reality, different scales grow at different 
rates, due to non-linearities. We show that angular surveys are sensitive to a weighted 
average of the power spectrum over a distribution of redshifts, where the weight func- 
tion varies with wavenumber. We compute this weight function and show that it is 
fairly sharply peaked at z, which is a function of k. As long as the extracted power 
spectrum is understood to be P(k,z), the error introduced by non-linear scaling is 
quantifiable and small. We study these effects in the context of the APM and SDSS 
photometric surveys. In general the weight matrix is peaked at larger z and is broader 
for deeper surveys, leading to larger (but still quantifiable) errors due to non- linear 
scaling. The tools introduced here - in particular the weight function and effective 
redshift z - can also be profitably applied to plan surveys to study the evolution of 
the power spectrum. 



1 INTRODUCTION 



Even as redshift surveys which allow us to obtain three di- 
mensional maps of the sky advance, photometric surveys still 
maintain their usefulness for cosmology. The fundamental 
advatange of the two dimensional surveys is that they can 
measure the positions of many more galaxies than can red- 
shift surveys. This is often enough to offset the loss of radial 
information, especially if one is interested only in some sim- 
ple statistics characterizing the underlying density field. 

In order to make sense of the angular information, one 
needs to understand how structure is sampled along the line 
of sight. For example, a deep survey picks up information 
about structure at much earlier times and much larger scales 
than a shallow one. Most of this information is encoded in 
the kernel which is given by Limber's Equation if the selec- 
tion function is known. Recapturing the three-dimensional 
power spectrum from the angular correlation function then 
involves inverting the kernel. This inversion is not com- 
pletely straightforward, but several different techniques have 
been used (Baugh & Efstathiou 1993,1994; Gaztanaga & 
Baugh 1998; Dodelson & Gaztanaga 1999; and Eisenstein & 
Zaldarriaga 1999), and they all seem to agree fairly well. 

One aspect of Limber's equation which is typically given 
short shrift is the question of how the power spectrum (or 
its Fourier transform, the correlation function) evolves with 
time. Some assumption is needed in order to generate the 
kernel; typically it has been assumed that the power spec- 
trum scales as (1 + z)~^ where z is the redshift and (3 = 2 
corresponding to linear evolution in a flat universe is the 



standard choice. It is important to note that making a choice 
is crucial to the success of the inversion process. An assump- 
tion about the time dependence of the power spectrum al- 
lows the kernel to be written as a function of wavenumber 
k only; undoing the integral over k for many angles is then 
possible. If no assumption was made, the integral would be 
over both redshift and k. It would be much more difficult, 
if not impossible, to undo this two-dimensional integration 
and get out P(k,z). 

Since we are forced into an assumption about how the 
power spectrum evolves with time, we need to ask how much 
this assumption affects the results. Here we examine this 
question. Section 2 briefly reviews the standard derivation 
of Limber's Equation. Section 3 introduces a tool to analyze 
the effectiveness of the assumption that P oc (l + z) - ' 3 . This 
is the recent work which allows one to generates a full non- 
linear P(k, z) from a given linear power spectrum. Armed 
with this tool, we then show two power spectra, both of 
which give the same angular correlation function. One has 
the simple (l + z)~ l3 scaling, while the other has more realis- 
tic scaling accounting for non-linearities. Although these two 
3D power spectra give the same angular correlation func- 
tion, they are much different today. We illustrate this for an 
APM-like survey (Maddux et al. 1990) . 

Section 4 isolates the reason for the difference between 
the two spectra. Essentially, any survey is actually a measure 
of the power spectrum over a range of redshifts centered 
at z, an easily computable function of wavenumber. If one 
insists on interpreting the results as measures of the power 
spectrum today, different scalings from z to z = lead to 
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different P(k,z = 0). However, if one interprets the results 
as a measurement of P(k,z), the scaling scheme one uses 
is irrelevant. Another way of saying this is to emphasize 
that the measurement is of P(k,z); the weighting function 
is fairly compact and so is often insensitive to the behavior 
of P for z much different than z. 

Finally section 5 computes the error in the estimate 
of P(k, z) introduced by assuming linear scaling through z. 
This error is largest on small scales where non-linear evolu- 
tion sets in earliest. For 3D wave numbers k 1 h Mpc~ , 
the error is quite small for APM, and larger, but still less 
than the statistical errors for a wider, deeper survey such as 
the Sloan Digital Sky Survey (SDSS)Q . 



2 THE KERNEL AND THE STANDARD 
ASSUMPTION 

We begin with the discretized, relativistic version of Lim- 
ber's equation, 



= V/c iaa Pao 



(i) 



where Wi is the angular correlation function in a bin cen- 
tered at Qi\ P aa is the power spectrum in a bin centered at 
wavenumber k a and redshift z a \ and K, is the kernel which 
depends on all three variables. The kernel is 

(A*)(A*) k a x(z a ) 4 



2tt (1 + z a ) e F(z a )H(z a ) 
x^ 2 (z a )Jo(k a x(z a ) sin 9i), 



(2) 



where x is the comoving distance out to redshift z; F de- 
pends on the cosmological model, equal to one in a flat, 
matter dominated universe; H is the Hubble rate as a func- 
tion of redshift which is also model dependent; and t[> is the 
selection function. It has been normalized so that 



dxx 2 a 3 (x)4>(x) = 1. 



(3) 



If one assumes a linearly evolving power spectrum, then 



P(k,z) =P(k,z = 0)(l + z)~ 



(4) 



and P(k, 0) can be moved out of the sum over redshifts. In 
this case, the kernel simplifies and we are left with: 



K P 



(5) 



where P a now denotes the power spectrum today and the 
new kernel is independent of redshift: 



K ia 



(Afc)(Az) 
2^ 



x(z a ) 4 



*F{z a )H{z 



xip 2 (z a )J (k a x(z a )smOi). (6) 

Equation (^) is then inverted to extract the three dimen- 
sional power spectrum P a . 

To reiterate, the separation in eq. ^ is certainly not 
correct, since the power on small scales evolves differently 
over time than the power at large scales. We now turn to 
more realistic scaling. 
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Figure 1. Present day values of three dimensional power spectra 
which might be inferred from angular data. The solid line (Pi ) is 
the present value of a spectrum which (i) evolves linearly and (ii) 
fits the angular data. The dashed and dotted curves (P2 , APM 
and P2, SDSS) are the present value of a spectrum which evolves 
non-linearly (i.e. realistically) and fits the agnular data. The dot- 
dashed curve (P3) is the linear spectrum which gives rise to the 
non-linear spectrum P2. It should be used to fit cosmological pa- 
rameters. 



3 NON-LINEAR SCALING 

To arrive at a non-linear power spectrum from a linear one, 
we use the treatment described by Peacock & Dodds (1996). 
In their paper they work in terms of the dimensionless power 
spectrum A 2 , where 



A 2 (k) = — P(k). 



(7) 



They introduce the non- linear wave- number, fcjvt, a func- 
tion of the linear wave-number, fci and the non- linear power 
In particular, 



spectrum A 2 NL 



kNL 



and 



(1 + A 2 NL (k NL )) 3 k L 



A 2 NL (k NL ) = f NL (A 2 L (k L )), 



(8) 



(9) 



where the function /nl is given in Peacock and Dodds. 

Armed with this transformation, we can invert it using 
the Newton- Raphson method to determine the linear power 
spectrum that would give rise to a given non-linear one. 
This linear power spectrum can be evolved trivially from 
early times until today, at each step of the way using eq. Q 
to form the corresponding non-linear power spectrum. This 
gives a much more realistic P(k, z), which can then be used 
to compute the angular correlation function. 

Figure [l] shows three power spectra which might con- 
ceivably be extracted from the APM survey. The three lines 
correspond to three different power spectra, all shown at 
z = 0: 
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Figure 2. Measured w used to fit the three-dimensional power 
spectrum and w calculated from the fit power spectrum using the 
non-linear method described herein. 



• Pi is the power spectrum one gets from the inversion 
assuming linear scaling. 

• Pi has the more realistic scaling using the formalism of 
Peacock and Dodds, but leads to a very similar w(8) (see 
figure 2). We will discuss in the next section how we arrived 
at this power spectrum. 

• Finally P3 is the linear spectrum associated with the 
non- linear spectrum P%. That is, if the universe started with 
P3 at early times (scaled back by (1 + z)~ 2 ), the non-linear 
power today would be Pi. 

The first of these, Pi, is what emerges from a blind in- 
version assuming P oc (1 + z)~ 2 . This is what is usually 
used to compare with theories. The second, Pi, is a much 
more accurate extraction of the power spectrum today since 
we accounted for the non-linear evolution. It clearly differs 
from Pi at k <; 0.5, so any attempt to use Pi to constrain 
cosmological parameters will necessarily be inaccurate. The 
only problem with P2 is that we have not yet explained how 
we got it. We'll do this in the next section. To do accurate 
parameter estimation, one does a best fit to the data al- 
lowing for several free parameters and evaluating the power 
spectrum at any point in parameter space with for example 
the BBKS (Bardeen et al. 1986) form or the output from 
CMBFAST (Seljak & Zaldarriaga 1996). To do this com- 
parison properly, the power spectrum P3 would need to be 
used, for the codes compute the linear power spectra and P3 
is the linear spectra corresponding to Pi. It has been com- 
mon practice to neglect non-linear effects and simply use Pi 
to fit for cosmological parameters, neglecting information 
on scales larger than some fc max , typically chosen to be in 
the range 0.1 — 0.2h Mpc -1 . Since Pi (the incorrect linear 
spectrum) differs from P3 (the correct linear spectrum) at 
wavenumbers k even smaller than O.l/i Mpc -1 , it would be 
clearly be much better to find a systematic way of obtaining 
P3 or equivalently its non- linear counterpart, P2. 




k = 0.01 



k = 0.1 



k=1.0 



Figure 3. The weight function W(k,z) for k = 1.0, 0.1, 0.01ft 
Mpc -1 for both the APM survey and SDSS, which goes wider 
and deeper. This shows which values of z contribute the most 
weight in the determination of the angular correlation function. 
The 3D power spectrum inferred therefore is P(k, z) at the peak 
of the wave function. 



4 THE WEIGHT FUNCTION 

We are almost ready to divulge the secret of how we got the 
spectrum Pi, which we claim is a much better estimate of 
the power spectrum today than is Pi. First, though, let's 
try to understand why the spectra Pi and Pi above differ. 
If the measurement was only of the power spectrum today, 
then it wouldn't matter how the spectrum evolved with time 
beforehand: to fit the angular correlation function, Pi would 
have to be equal to Pi today. The difference between the two 
today, as reflected in figures 1 and 2, then must be due to 
the fact that the measurements are a weighted average of 
the power spectrum over time. Different surveys will carry 
with them different weight functions. Indeed, even different 
wavenumbers in the same survey will have different weight 
functions. It is clearly very important to understand and be 
able to compute the weight function. 

The weight function can be computed by first forming 
a x 2 from the observed un's and the theoretical Wi(P): 

X 2 = (Wi - Wi(P)) C-/ (w 3 - Wj (P)) (10) 

where C is the covariance matrix for Wi . If we want to figure 
out how much weight the power spectrum at redshift z a and 
wavenumber k a contributes to the x 2 , we need only compute 
the second derivative of the x 2 with respect to P aa . This 
gives the curvature, or the weight function: 



W{k a ,z a 



15V 



2 dPL 



ICi, 



*Cy ICjaa- 



(11) 
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Plugging in from eq. (H), we see that the weight function is 



W{k,z) = 



k(Ak)(Az) 



2tt HF(1 + z) 
x C^j 1 Jo{kxsm6i)Jo(kxs'mdi). 



(12) 



for two sur- 
For the for- 



The weight function is plotted in figure 3 
veys, APM and the SDSS photometric survey T 
mer, which is shallower, the weight function is peaked at 
z ~ 0.08 and is fairly narrow. The weight function for SDSS 
is also shown assuming galaxies can be extracted down to 
22nd magnitude. As expected it peaks at higher redshift and 
is broader. For fixed k we define z to be the redshift at which 
W(k, z) peaks. 

The weight functions shown in figure ^| may be some- 
what suprising to those with knowledge of the surveys. The 
median redshifts of these surveys are larger than might be 
expected from consideration of figure ^. The discrepency can 
be attributed to the fact that, for a given k, quite a bit of the 
weight for the measurement of the power spectrum comes 
from large angles (and therefore reshifts much smaller than 
the median redshift). 

The weight function gives us a very clear way to think 
of the power spectra extracted from the inversion. Recall 
that the inversion techniques assume linear scaling. Since 
the weight function tells us that a given k— mode is mostly a 
measure of P(k, z), we should scale back the inverted power 
spectrum to z. Then, if we want the power spectrum today, 
we can scale forward with the non-linear formulae. Indeed, 
we can now reveal that this is how we arrived at Pi in the 
previous section. 

This suggests the following recipe for extracting a 
present day 3D power spectrum from angular data: 

• Assume linear scaling so that Eq. (1) reduces to the 
much more managable Eq. (5). 

• Invert to find Pi today. 

• Scale P\(k,z — 0) back linearly to Pi(k,z(k)) = 
Pi(k, z — 0)(1 + z(k))~ 2 . This scaled back spectrum is a 
good estimator for the non-linear spectrum at z. 

• Scale the spectrum obtained in the previous step non- 
linearly to its present value, Pi. 

• Find the underlying linear spectrum corresponding to 
Pi and Pz, call it Pz- This can then be compared to linear 
models to extract parameters. 



5 ERROR ON P(K, z) 

There is one final loose end to tie up. The above prescription 
would be exact if the weight function was a delta-function, 
infinitely sharp at z. Its finite width allows for the possibility 
that the scaling assumed around z affects the measurement 
of P(k, z). There are several ways we can test this possibil- 
ity. The first is to look at the resultant angular correlation 



t To do this, we have had to assume something about the co- 
variance matrix. We restricted ourselves to angular scales greater 
than half a degree and assumed the covariance matrix was due 
solely to cosmic variance. We computed this matrix assuming 
Gaussian statistics. 
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a measure of the error induced by assuming linear scaling to the 
statistical error. 



functions from the two spectra. Figure 2 shows these and the 
difference between the two for APM. It is encouraging that 
the difference is so small. This suggests that the evolution 
of the power spectrum through z is not very important; the 
measurements are simply of P(k, z). 

To test this further and to assign error bars due to the 
assumed scaling, we can define 



(13) 



Since ^ W(k,z a ) is the weight of the measurement, it is 
the inverse of the square of the error on the measurement 
of P(k). Therefore, Asc can be thought of as the ratio of 
the error due to the linear scaling assumption used in the 
inversion process to the overall error in the power spectrum. 
If Asc is much less than one, then we need not worry about 
the scaling assumption. Figure ^ shows that, for surveys 
with the depth of APM or even SDSS, Asc is indeed quite 
small for k Ih Mpc -1 . The broader weight function of 
SDSS leads to a larger scaling error. One might argue that, 
even for k as low as 0.3h Mpc -1 , the statistical error on 
the power spectrum is an underestimate. Taking As C to be 
0.2 there leads suggests (assuming errors add in quadrature) 
that the linear evolution assumption increases the errors by 
a factor of y/1 + 0.2, about ten percent. And of course, at 
higher wavenumbers the error bars get even larger. 



6 CONCLUSIONS 

The inversion of the angular correlation function gives a 
measure of the three-dimensional power spectrum at red- 
shift z, where z is the place where W(k, z) in eq. (^) peaks. 
A simple way to obtain an estimate for P(k, z) is to assume 
linear scaling of the power spectrum, invert the kernel, and 
scale back the power spectrum to z. An estimate in the error 
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incurred by this procedure is given by eq. (|l3|). For current, 
and even future surveys, the error is small for k < 0.5h 
Mpc - (but might be much more significant for deep sur- 
veys). On larger scales, the error is larger. This is not nec- 
essarily a bad thing: it is an indication that the survey is 
sensitive to the evolution of the power spectrum. In fact, 
the tools developed here could be applied to help plan sur- 
veys or devise optimal strategies for breaking a survey into 
subsets. One could compute the weight function for a given 
subset of data (e.g. a given magnitude slice) and choose a 
different subset whose weight function peaks far away in 
redshift. The width of the weight function in SDSS suggests 
that this may already be possible. 

We have not dealt at all with the possibility of using 
photometric redshifts to learn more about the evolution of 
the power spectrum. And we have completely ignored the 
issue of how the galaxies are biased with respect to the mat- 
ter. There has been much activity in both of these fields over 
the past several years, which should help extract even more 
useful information from angular surveys. 
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